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At present, there is a worldwide effort to use cold atoms to simulate strongly correlated quantum 
many-body systems. It is hoped that these "simulations" will provide solutions to many unsolved 
problems. However, the relevant energy scales in most of these experiments are so small that one 
has to go to entropy regimes far below those achievable today. Here, we present a general scheme to 
extract entropy directly from the region of interest. The late stage of this process is equivalent to a 
continuous "evaporation" , and is able to combat intrinsic heating of the system. For illustration, we 
show how to cool a weak coupling BCS superfluid (with T c ~ 10 _10 K) to 10 _11 K with this simple 
procedure, with entropy per particle as low as 5 x lO -4 /^ in the superfluid region. 



Although quantum simulations using cold atoms are 
considered as means to find solutions to unsolved many- 
body problems, they are in fact an exceedingly ambitious 
undertaking in ultra-low temperature physics. For exam- 
ple, the antiferromagnetic phase of the 3D fermion Hub- 
bard model is characterized by the virtual hopping scale 
t 2 jJJ ', which can be 10 _11 K or lowerpQ. In terms of en- 
tropy per particle, it is estimated that one has to lower 
the current experimentally achievable value by about a 
factor of 50 to reach the antiferromagnetic phase not too 
far below T c [2]. An even larger factor is required to reach 
the low temperature regime. In addition, intrinsic heat- 
ing rate due to spontaneous emission can lead to entropy 
production as high as lfc^ per second [3 , which will de- 
stroy many states of interest over the duration of the ex- 
periment. It is therefore crucial to invent cooling schemes 
powerful enough to overcompensate the intrinsic heating. 

The 3D fermion Hubbard model, however, is by no 
means the most demanding for temperature and entropy 
reduction. There is a host of novel cold atom states 
characterized by very small energy scales. For example, 
high spin fermions such as 40 K have ten Fermi surfaces 
and can have very remarkable pairing states in optical 
traps [4 j. Unfortunately, the interactions are so weak 
that the superfluid transition temperature can be as low 
as 10 _11 K for systems with 10 5 particles, making the re- 
alization of such novel states very difficult. Other systems 
such as large spin Bose gases in optical lattices, mixtures 
of quantum gases, and low dimensional dipolar Fermi 
gases [5 are characterized by similar or even smaller en- 
ergy scales. Thus, the success of realizing many novel 
states rests on the ability to reach these unprecedented 
ultra-low temperature and entropy regimes. 

Previously, we proposed a scheme to expel (or 
"squeeze" out) the entropy of a Fermi gas in an optical 
lattice to the surface of the cloud by turning its interior 
into a band insulator through compression [6]. The sur- 
face entropy can then be removed by various means. For 
example, the surface entropy can be transferred into a 



surrounding Bose-Einstein condensate (BEC) which can 
then be evaporated away [6]. It is shown that in this way, 
the entropy per particle can be lowered by a factor of 
50. Other variations of using band insulators to push 
out the entropy have been considered. [HE]. The effec- 
tive transfer of entropy between quantum gases has also 
been demonstrated recently [9]. It was also pointed out 
in Ref. [6] that this principle of entropy expulsion is not 
restricted to band insulators, but applicable to all system 
with an incompressible gapful phase. 

The "expulsion method" [6] and its variations [7J [8] , 
however, have many limitations. First of all, many sys- 
tems of interests do not have a band insulator or an in- 
compressible phase. This include all quantum gases in 
single traps without an optical lattice, as well as spinor 
gases and mixtures in optical lattices. Secondly, it in- 
volves too many steps, requiring first the creation of a 
band insulator, then the removal of entropy at the sur- 
face in a way that will not lead to entropy re-generation, 
and then the transformation of the band insulator into 
the state of interest. Thirdly, its cooling power, though 
considerable, is still limited. It is not clear whether it 
can overcome the problem of intrinsic heating. While in 
principle one can repeat this process, the number of steps 
required makes a repeated running impractical. 

In this paper, we present a cooling scheme free of the 
above limitations. The scheme is applicable to all quan- 
tum gases including mixtures and gases with internal de- 
grees of freedom. It is designed to extract entropy di- 
rectly from the region of interest in a continuous way, 
and hence is capable of combating extrinsic and intrinsic 
heating in a continuous manner. This cooling scheme will 
remain operative as long as the system can equilibrate 
during the operation. However, as temperature drops, 
the time for equilibration increases and will eventually 
place a limit on the cooling power of this scheme. Here, 
we focus on ways to reach maximum cooling power within 
adiabatic processes. How to speed up the equilibration 
processes will be considered elsewhere. 
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FIG. 1: Schematic of the entropy extraction scheme: The sys- 
tem of interest is confined in A. The trapping potential V(x) 
is represented by the solid blue curve. Chemical potential \i 
is represented the red horizontal line. As the surrounding re- 
gion B expands, entropy and particles flow out from A to B. 
/jL is kept close to the top of the dimple, and the difference 
jl = /jl — V(0) is maintained roughly constant at the late stage 
of expansion so that T/p, decreases as expansion proceeds. 



A. Direct entropy extraction: The general prin- 
ciple for our scheme is shown in Figures (1A) to (1C). 
The system of interest is confined in region A. It is sur- 
rounded by an "extraction region", denoted as B. An 
example of the trapping potential V(x) that gives rise 
to this situation is shown in Figure (IB). It consists of a 
box-like potential (region B) with an attractive "dimple" 
at the center (region A). Region B is further divided 
into regions B\ and B2, where \i > V(x) and /i < V(x) 
respectively. If we expand the volume of B adiabati- 
cally while keeping that the chemical potential fi close 
to the top of the dimple (Fig.(lC)), both energetic par- 
ticles and entropy will be "sucked out" from A to £>, 
causing the temperature T to drop. More importantly, 
if the flow entropy flow out of A is faster than the parti- 
cle flow, then the entropy per particle in A (denoted as 
(S / N) dimple) a l so decrease. Next, we note that the 
entropy per particle in A is proportional to T/jl, where 
jl = \i — V(0) is the chemical potential measured from 
the bottom of the dimple, (see example below). One can 
therefore reduce (S/N) dimple as B expands by keeping jl 
at a reasonably high value. In addition, we note that 
in the low-fugacity region 62 > the entropy per particle is 
s/n = fc# [5/2 — ///(fceT)], where s and n are entropy 
and particle density, respectively. In this region, we have 
n = e ^/^ T )/A 3 < 1, where A = h/^27rMk B T is the 
thermal wavelength. Since in B2, /i(x) < and /i(x) in- 
creases as |x| increases, the entropy per particle s(r)/n(r) 
becomes arbitrarily large as one moves away from the 
center, making B2 a highly efficiently region to absorb 



FIG. 2: Trapping potential for Fermi gas with weak attractive 
interaction in the z — plane. The blue curve is the total 
potential which consists of a dimple and the harmonic poten- 
tial. The horizontal red line denotes the chemical potential. 
The inset shows the potential on a larger scale. 



entropy. Our discussion of a box-like potential is for il- 
lustration purposes only. In practice, one simply uses a 
harmonic trap. The expansion of B then corresponds to 
reducing the trap frequency. 

We should mention that the dimple trap has been used 
before by Ketterle's group to produce a Bose-Einstein 
condensate (BEC) [TT . In that case, the dimple is turned 
on without changing the harmonic trap. This actually 
raises the temperature [TT]. Nevertheless, BEC still oc- 
curs because Tbec rises even faster due to the rapid in- 
crease in phase space density. This scheme, however, 
does not lead to large reduction of entropy per particle 
within the dimple region. To achieve that, it is crucial to 
open up the harmonic trap while controlling the density 
at the center of the dimple and the value of the chemical 
potential as we show below. 

B. Illustration of the Cooling Scheme: To il- 
lustrate the cooling power of this scheme, we consider 
a two-component Fermi gas with weak attractive inter- 
actions. At present, superfluidity in Fermi gases can 
only be achieved in the strongly interacting limit with 
|fcpa s | ^> 1, where kp is the Fermi wave vector and a s is 
the s- wave scattering length. For a weakly attractive gas, 
\kpa\ <C 1, the superfluid transition temperature T c is too 
low to be reached by the cooling methods today. Here, we 
show that such ultra-low temperatures are achievable by 
our simple scheme. We begin with a (two-component) 
Fermi gas of N = 1 x 10 6 6 Li atoms with scattering 
length a s = — 0.25 /hp- The corresponding T c in the bulk 
is T c = 0ME F e- 7r ^ 2kFas \ [10 , where E F = h 2 k 2 /2M is 
the Fermi energy, (for kp = 2.2/im -1 , T c = 0.24nK. ) 
Our trapping potential is V(x) = V^(x) + Vd( x )> where 
Vh and Vd are the potentials of a harmonic trap and a 
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calculate these quantities and the temperature within the 
local density approximation as V(r) is varied. The prop- 
erties of the system in the final potential shown in Figure 
2 are given by Figures 3 to 5, as well as in Table I. 
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FIG. 3: Density and energy gap of the Fermi gas in the po- 
tential shown in Figure 2. The parameters of this system are 
given in Table 1. 



dimple trap respectively, 



^( x ) = Y^ ±r 



y d (x) = V (l - e~ r ' D ). 

(1) 

Here, r = (x,y), x = (x,y,z); co_\_ and uo z are the fre- 
quencies of the harmonic trap; V Q and D are the depth 
and the width of the dimple. We have chosen a dimple 
with no z- variations to mimic a red de-tuned laser along 
z. (The specific form of the dimple, however, is not im- 
portant for our discussions). We have added a constant 
to the dimple potential so that the minimum value of the 
total potential is zero. 

We consider an initial state where the Fermi gas is 
confined in a harmonic trap without the dimple, 

V®(r) = ^M^ 2 (r 2 + z 2 ), V d = 0; (2) 

with entropy per particle (S/N)^ = 0.45, which is close 
to the lowest value achievable in a single trap today. For 
/2tt = 23.7Hz, this entropy corresponds to an 
initial temperature Ti = 9AnK, and a Fermi momentum 
at the center of the trap kp = 2.2jim~ 1 . Next, we turn 
on the dimple and open up the harmonic trap adiabat- 
ically. The depth and width of the dimple are adjusted 
so that density (and hence kp) at the center of the dim- 
ple remain constant during the entire process. The final 
configuration of the trapping potential is shown in 
Figure 2 at two different length scales. The potential 
depth V is measured in units of the Fermi energy Ep 
at the center of the dimple, which is fixed to its initial 
value Ep = Ti 2 k 2 F /2M. The frequency of the dimple trap 
near its bottom is 361 Hz. (See also Table 1.) The size of 
the dimple is D = 10.6/im. The horizontal straight line 
indicates the value of the chemical potential \i. 

Using weak-coupling BCS theory, one can obtain the 
number density n, entropy density s, and energy gap A 
as functions of \i and T. (See Appendix.) We can then 



Figure 3 shows the density profile and the gap. For 
these trap parameters, the temperature has been reduced 
to one-fifth T c , and pair condensation has taken place 
inside the dimple. Figure 4 shows the entropy density 
s(r) and entropy per particle s(r)/n(r) as a function 
of r in the z = plane at two different scales. The 
high entropy density accumulated at the edge of the 
dimple is due to the rapid drop of density there. The 
entropy per particle in the dimple region remains very 
low. It is more illuminating to look at the total particle 
number N(R) and total entropy S(R) contained within 
a region of radius R, S(R) = (27rrdr) / dzs(r, z), 

N(R) = J^(27rrdr) J dzn(r, z); which are plotted in Fig- 
ure 5. We find that the dimple (R < D) contains about 
0.15 of total number of particles, whereas essentially all 
the entropy resides outside the dimple. Counting the 
total entropy and particle number in the dimple region 
(consisting of both normal and superfluid components) , 
we find that their ratio is (S/N)^ = 0.0014, reduced 
from the original value by a factor of 320. Counting 
only the superfluid region within the dimple, the ratio 
(S/N)faln sf 1S 900 lower than the initial value, (see Ta- 
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FIG. 4: Entropy density s(r) and local entropy per particle 
s(r)/n(r) for the Fermi gas shown in Figure 3 are represented 
by the red and blue curves. The inset is the same figure in 
much larger scale, with the same labeling of the axes. 
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FIG. 5: Total particle number and total entropy within a 
distance R from the center. 

ble 1). We stress that this reduction is by no means the 
limit. More cooling can be obtained by further opening 
the harmonic trap and adjusting the dimple potential. 

C. Equivalence to continuous evaporation: As 
the potential outside the dimple becomes flatter and flat- 
ter, or more generally, as the volume of the region B in 
Figure 1 continues to grow, the system inside the dimple 
(region A in Figure 1 ) can not tell whether the potential 
outside is expanding but eventually bounded, or a flat 
potential that turns downward at large distances, allow- 
ing particles to leak out to infinity. In other words, as the 
trap opens up adiabatically, as far as the system in A is 
concerned, it amounts to a process of continuous evapora- 
tion. From this viewpoint, one can see another advantage 
of this scheme: the continuous removal of entropy in A 
allows one to combat extrinsic or intrinsic heating. By 
continuous adjusting the depth V Q and width D of the 
dimple, one can keep the chemical potential close to the 
top of the dimple, therefore continuously "evaporating" 
away the hot particles inside. 

Our scheme differs, however, from the usual evapo- 
rative cooling in a significant way. In the evaporation 
scheme used by current experiments, one lowers the po- 
tential near the edge of the cloud to let the hottest par- 
ticles escape [12 . While temperature is lowered in this 
process, the density at the center of the trap (and hence 
Ep) is also decreased. As a result, once T falls below Ep, 
there is no significant reduction in the ratio of T/Ep [12] , 
and hence no significant reduction on entropy per par- 
ticle in the dimple, since (S/N) ~ T/Ep for a Fermi 
gas. In our case, the density at the center of the dim- 
ple is kept fixed during the evaporation process. The ratio 
T/Ep drops rapidly as the cooling process proceeds. This 
simple step makes an enormous difference in the cooling 
power of the evaporation process. 

Our discussions assume that the process is adiabatic. 
As mentioned in the opening, as temperature drops, it 
will take the system in the dimple longer and longer to 
equilibrate. Eventually, it will reach a point where adia- 
baticity can not be maintained. The equilibration time, 



which is system-dependent, is the intrinsic limit of the 
cooling power of this scheme. The speeding up of relax- 
ation processes will be studied elsewhere. 

Finally, we point out that the principle discussed in A 
applies to all quantum gas systems. For lattice quantum 
gases, Eq. will be a potential in addition to the lattice 
potential. The crucial steps remains the same: relaxing 
the harmonic trap adiabatically, adjusting the dimple so 
that fx remains close to the top, and jl is kept at a value 
necessary to realize the state of interest. Despite the 
simplicity of these steps, they make big differences from 
previous studies [TTJ [12]. The physics contained in these 
steps can work magic, enabling entropy per particle to be 
lowered from the current values by orders of magnitude. 

Appendix: Within BCS theory, number den- 
sity n and entropy density s are given by n = 

(* - = £k-M, £ k =V£k + A2 > £ k = h 2 k 2 /(2M), 

E k [(l " /(£k))ln(l - f(E k )) + f(E k )lnf(E k )}, 
where fi is the volume of the system, f(x) = l/(e x + 1) 
is the Fermi function, and the gap A is determined by 

M 1 ^ f tanhE k /(2k B T) 1_\ 

4h 2 ira 2^\ E k ej- W 

Within the local density approximation, the total particle 
number N and total entropy S are given by TV = J n(ja — 
F(x,T)) and S = J s(fi - V(x,T)). For given values of 
N and 5, one can then find \i and T, which then gives 
the results in Figs. 3-5. 
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